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Abstract 


Transverse momentum spectra of charged particles are measured by the CMS ex¬ 
periment at the CERN LHC in pPb collisions at ~ 5.02 TeV, in the range 

0.4 < pj < 120GeV/c and pseudorapidity |?/cm| < 1.8 in the proton-nucleon center- 
of-mass frame. For pj < lOGeV/c, the charged-particle production is asymmetric 
about rjcu = 0, with smaller yield observed in the direction of the proton beam, qual¬ 
itatively consistent with expectations from shadowing in nuclear parton distribution 
functions (nPDF). A pp reference spectrum at ^Js = 5.02 TeV is obtained by inter¬ 
polation from previous measurements at higher and lower center-of-mass energies. 

The Pt distribution measured in pPb collisions shows an enhancement of charged 
particles with pj > 20 GeV/c compared to expectations from the pp reference. The 
enhancement is larger than predicted by perturbative quantum chromod}mamics cal¬ 
culations that include antishadowing modifications of nPDFs. 
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1 Introduction 


The central goal of the heavy ion experimental program at ultra-relativistic energies is to cre¬ 
ate a system of deconfrned quarks and gluons, known as a quark-gluon plasma (QGP), and 
to study its properties as it cools down and transitions into a hadron gas. A key tool in the 
studies of the QGP is the phenomenon of jet quenching ll], in which the partons produced in 
hard scatterings lose energy through gluon radiation and elastic scattering in the hot and dense 
partonic medium |2J. Since high transverse momentum quarks and gluons fragment into jets 
of hadrons, one of the observable consequences of parton energy loss is the suppression of the 
yield of high-px particles in comparison to their production in proton-proton (pp) collisions. 
This suppression, studied as a function of the px and pseudorapidity (p) of the produced par¬ 
ticle, is usually quantified in terms of the nuclear modification factor, defined as 


Rab(pt/ p) 


1 d^N^^/dpxdp 
(Tab) dVPP/dpxdp ' 


( 1 ) 


where is the particle yield in a collision between nuclear species A and B, crPP is the cor¬ 
responding cross section in pp collisions, and (Tab) is the average nuclear overlap function 0| 
in the AB collision (in the case of proton-nucleus collisions, the quantity (Tab) = (TpA) is 
called average nuclear thickness function). If nuclear collisions behave as incoherent super¬ 
positions of nucleon-nucleon collisions, a ratio of unity is expected. Departures from unity 
are indicative of final-state effects such as parton energy loss, and/or initial-state effects such 
as modifications of the nuclear parton distribution functions (nPDF) [^. The nPDFs are con¬ 
strained by measurements in lepton-nucleus deep-inelastic scattering (DIS) and Drell-Yan (DY) 
production of dilepton pairs from qq annihilation in proton-nucleus collisions |5i|. In the small 
parton fractional momentum regime (x < 0.01), the nPDFs are found to be suppressed rela¬ 
tive to the proton PDFs, a phenomenon commonly referred to as "shadowing" |6|. At small 
X, where the parton distributions are described theoretically by non-linear evolution equa¬ 
tions in X, gluon saturation is predicted by the color glass condensate models 10-El- For the x 
regime 0.02 < x < 0.2, the nPDFs are enhanced ("antishadowrng") relative to the free-nucleon 
PDFs |5i 

To gain access to the properties of the QGP produced in heavy ion collisions it is necessary to 
separate the effects directly related to the hot partonic medium from those that are not, referred 
to as "cold nuclear matter" effects. In particular, nPDF effects are expected to play an important 
role in the interpretation of nuclear modification factors at the GERN LHG. Unfortunately, the 
existing nuclear DIS and DY measurements constrain only poorly the gluon distributions over 
much of the kinematic range of interest. High-px hadron production in proton-nucleus (or 
deuteron-nucleus) collisions provides a valuable reference for nucleus-nucleus collisions, as 
it probes initial-state nPDF modifications over a wide kinematic range and is expected to be 
largely free from the final-state effects that accompany QGP production IITOl . 

The measurements of the nuclear modification factors of neutral pions and charged hadrons in 
the most central gold-gold (AuAu) collisions at the Relativistic Heavy Ion Gollider (RHIG) lITT] - 
[m revealed a large suppression at high pj, reaching an Raa as low as 0.2. In contrast, no such 
suppression was found at mid-rapidity in deuteron-gold collisions at the same energy IITST 
[TSl . These findings established parton energy loss, rather than initial-state effects dll, as the 
mechanism responsible for the modifications observed in AuAu collisions. 

At the LHG, the charged-particle suppression in lead-lead (PbPb) collisions persists at least up 
to a Px of lOOGeV/c ||20H^. In proton-lead (pPb) collisions, the ALIGE Gollaboration reported 
no significant deviations from unity in the charged-particle Rppb up to px ~ 50GeV/c Il22l . The 
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2 Data selection and analysis 


analysis presented here used data from CMS to extend the measurement of the charged-particle 
RpPb out to Pt ~ 120GeV/c, with the aim of evaluating initial-state effects over a kinematic 
range similar to that explored through measurements in PbPb collisions Il20l . 

Proton-nucleus collisions have already been used to assess the impact of cold nuclear matter 
on jet production at the LHC. The transverse momentum balance, azimuthal angle correlations, 
and pseudorapidity distributions of dijets have been measured as a function of the event activ¬ 
ity, and no significant indication of jet quenching was found [231. When normalized to unity, 
the minimum-bias dijet pseudorapidity distributions are found to be consistent with next-to- 
leading-order (NLO) perturbative quantum chromodynamic (pQCD) calculations only if nPDF 
modifications are included [241 . Similarly, inclusive jet Rppb measurements are also found to be 
consistent with NLO pQCD predictions that include nPDF modifications Il25ll26[ . The measure¬ 
ment of the charged-particle spectra presented in this paper provides a comparison to theory 
that is sensitive to smaller x values than those accessible in the jet measurements. However, 
it should be noted that the charged-particle Rppb is dependent upon non-perturbative hadron- 
ization effects, some of which, such as gluon fragmentation into charged hadrons, are poorly 
constrained at the LHC energies 127]. 

The Pj distributions of inclusive charged particles in pPb collisions at a nucleon-nucleon center- 
of-mass energy of 5.02 TeV are presented in the range of 0.4 < pj < 120GeV/c. The mea¬ 
surement is performed in several pseudorapidity intervals over |?/cm| < 1-8. Here //cm is the 
pseudorapidity in the proton-nucleon center-of-mass frame. The nuclear modification factor is 
studied at mid-rapidity, |//cm| < 1/ and the forward-backward asymmetry of the yields, Yasym/ 
defined as 

M), . /r;d//cMd2Nch(pT)/d//cMdpT 

Y^sym(PT) = —T-;-^^, (2) 

d//cM d2Nch(pT)/d//cM dpx 

is presented for three pseudorapidity intervals, where a and b are positive numbers, and N^h is 
the yield of charged particles. 

Due to their wide kinematic coverage, the measurements are expected to provide information 
about the nPDFs in both the shadowing and antishadowing regions. In particular, the effects 
of shadowing are expected to be more prominent at forward pseudorapidities (in the proton¬ 
going direction), where smaller x fractions in the nucleus are accessed. 

In the absence of other competing effects, shadowing in the Pb nPDFs would result in values 
of Yasym > 1 at low Pj (i.e., small x). The effects of antishadowing can be probed with the Rppb 
measurement at larger px values of 30 < px < 100 GeV/c that correspond approximately to 
0.02 < X < 0.2. Antishadowing in the nPDFs may increase the yield of charged particles in pPb 
collisions compared with expectations from the yield in pp collisions. 


2 Data selection and analysis 

2.1 Experimental setup 

A detailed description of the GMS detector can be found in Ref. [28|. The GMS experiment 
uses a right-handed coordinate system, with the origin at the nominal interaction point (IP) 
at the center of the detector, and the z axis along the beam direction. The silicon tracker, lo¬ 
cated within the 3.8 T magnetic field of the superconducting solenoid, is used to reconstruct 
charged-particle tracks. Gonsisting of 1440 silicon pixel detector modules and 15148 silicon 
strip detector modules, totaling about 10 million silicon strips and 60 million pixels, the silicon 
tracker measures the tracks of charged particles within the pseudorapidity range |//| < 2.5. 
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It provides an impact parameter resolution of ~ 15 jim and a px resolution of about 1.5% for 
particles with px of 100 GeV/c. An electromagnetic calorimeter (ECAL) and a hadron calori¬ 
meter (HCAL) are also located inside the solenoid. The ECAL consists of more than 75 000 
lead tungstate crystals, arranged in a quasi-projective geometry; the crystals are distributed in 
a barrel region (|?/| < 1.48) and in two endcaps that extend out to |p| ~ 3.0. The HCAL barrel 
and endcaps, hadron sampling calorimeters composed of brass and scintillator plates, have an 
acceptance of |j/| < 3.0. The hadron forward calorimeters (HE), consisting of iron with quartz 
fibers read out by photomultipliers, extend the calorimeter coverage out to |p| = 5.2, and are 
used in offline event selection. Beam Pick-up Timing for the experiments (BPTX) devices were 
used to trigger the detector readout. They are located around the beam pipe at a distance of 
175 m from the IP on either side, and are designed to provide precise information on the LHC 
bunch structure and timing of the incoming beams. The detailed Monte Carlo (MC) simulation 
of the CMS detector response is based on Geant4 Il29| . 

This measurement is based on a data sample corresponding to an integrated luminosity of 
35 nb^^, collected by the CMS experiment in pPb collisions during the 2013 LHC running pe¬ 
riod. The center-of-mass energy per nucleon pair was = 5.02 TeV, corresponding to per- 
nucleon beam energies of 4 TeV and 1.58 TeV for protons and lead nuclei, respectively. The data 
were taken with two beam configurations. Initially, the Pb nuclei traveled in the counterclock¬ 
wise direction, while in the second data-taking period, the beam directions were reversed. Both 
data sets, the second one being larger by approximately 50%, were analyzed independently, 
yielding compatible results. To combine data from the two beam configurations, results from 
the first data-taking period are reflected along the z-axis, so that in the combined analysis, the 
proton travels in the positive z and t] directions. In this convention, massless particles emitted 
at rjcM = 0 in the nucleon-nucleon center-of-mass frame will be detected at = 0.465 in the 
laboratory frame. A symmetric region about rjcu = 0 is used in the data analysis, resulting in a 
selected pseudorapidity range of \t]cu\ < 1-8. 

2.2 Event selection 

The CMS online event selection employs a hardware-based level-1 (LI) trigger and a software- 
based high-level trigger (HLT). A minimum-bias sample is selected first by the LI require¬ 
ment of a pPb bunch crossing at the IP (as measured by the BPTX), and an HLT requirement 
of at least one reconstructed track with px > 0.4GeV/c in the pixel tracker. Lor most of the 
5.02 TeV data collection, the minimum-bias trigger is significantly prescaled because of the 
high instantaneous LHC luminosity. To increase the px reach of the measurement, a set of 
more selective triggers is also used: additional LI requirements are imposed to select events 
that have at least one reconstructed calorimeter jet with an uncorrected transverse energy of 
Ej > 12 GeV, and Ej > 16 GeV. These event selections are complemented by additional HLT 
requirements that select events based on the presence of at least one track with px > 12 GeV/c 
(for LI Ej > 12 GeV), or with px > 20 or 30GeV/c (for LI Ej > 16 GeV) reconstructed in the 
pixel and strip tracker. 

The above triggering strategy allows for the optimization of the data-taking rate while ade¬ 
quately sampling all px regions, including collecting all events containing very high-px tracks. 
The track trigger with a px threshold of 12 GeV/c records about 140 times more events with 
high-px tracks than the minimum-bias trigger, the track px > 20 GeV/c trigger enhances this 
with an additional factor of about 8, while the track px > 30 GeV/c trigger that is not prescaled, 
increases the number of events with a high-px track by yet another factor of about 2. 

In the offline analysis, additional requirements are applied. Events are accepted if they have 
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(i) at least one HF calorimeter tower on both the positive and negative sides of the HF with 
more than 3 GeV of total energy, (ii) at least one reconstructed collision vertex with two or more 
associated tracks, and (iii) a maximum distance of 15 cm along the beam axis between the vertex 
with the largest number of associated tracks and the nominal IP. Beam-related background is 
suppressed by rejecting events where less than 25% of all reconstructed tracks are of good 
quality ||30 |. 

An event-by-event weight factor accounts for correcting the measured charged-particle spectra 
in pPb collisions to a detector-independent class of collisions termed as "double-sided" (DS) 
events, which are very similar to those that pass the offline selection described above. A DS 
event is defined as a collision producing at least one particle in the pseudorapidity range —5 < 
'/lab < ~3 and another in the range 3 < J/iab < 5, each with proper lifetime t > 10^^^ s 
and energy E > 3 GeV ||3T]| . The performance of the minimum-bias and high-px single-track 
triggers, as well as the offline criteria in selecting DS events, is evaluated with simulations using 
the HIJING MG generator Il32ll . version 1.383 1331 . and the correction factors are computed as a 
function of the event multiplicity. An efficiency of 99% is obtained for the minimum-bias trigger 
and a negligible correction (i.e., 100% efficiency) for the high-px track-triggered events. The 
correction factor is also evaluated using an EPOS Il34l simulation and, based on the difference 
between both generators, a slightly px-dependent systematic uncertainty of 1% is assigned to 
the final spectra. 

During the pPb data taking period, about 3% of the recorded events contained more than one 
pPb collision. To reduce potential bias in the measurements arising from such "pileup", events 
with multiple reconstructed vertices are removed if the longitudinal distance between the ver¬ 
tices along the beamline is greater than a specific value that is related to the uncertainty of the 
vertex position. This value is also dependent on the number of tracks associated with each ver¬ 
tex and ranges from 0.2 cm for vertex pairs with at least 25 tracks associated with each vertex, 
to 3 cm for vertex pairs with only 3 tracks associated with the vertex having the fewest associ¬ 
ated tracks. Simulated HijiNG events are used to tune the pileup-rejection algorithm in order 
to reduce the number of erroneously eliminated single-collision events to a negligible fraction, 
and still maintain a high rejection efficiency for genuine pileup events. The pileup-rejection 
efficiency is found to be 92 ± 2%, which is confirmed by using a low bunch intensity control 
sample in data. 

To obtain inclusive particle spectra up to px ~ 120 GeV/c, data recorded with the minimum-bias 
and high-px track triggers must be combined appropriately. The corresponding weight factors 
are computed by counting the number of events that contain leading tracks (defined as the 
track with the highest px in the event) in the range of |?/iab | < 2.4 with px values in regions not 
affected by trigger thresholds, i.e., where the trigger efficiency of the higher-threshold trigger is 
constant relative to that of the lower-threshold trigger. The ratio of the number of such events 
in the two triggered sets of data are used as weight factors. For example, the region above 
which the track trigger with a px threshold of 12 GeV/c has constant efficiency is determined by 
comparing the px distribution of the leading tracks to that of the minimum-bias data. Similarly, 
the constant efficiency region of the 20 GeV/c track trigger is determined by comparison to the 
12 GeV/c track trigger, and the 30 GeV/c trigger to the 20 GeV/c trigger. The regions of constant 
efficiency for each trigger, as a function of leading charged-particle px, are shown in Fig. The 
12, 20, and 30GeV/c triggers have constant efficiencies above a leading charged-particle px of 
14, 22, and 32 GeV/c, respectively. The weight factors are then computed using the leading- 
track Px classes of 14 < px < 22GeV/c, 22 < px < 32GeV/c, and px > 32GeV/c for the 
three high-px triggers. The combined uncertainty in these normalizations is dominated by the 
matching of the 12 GeV/c track-triggered events to the minimum-bias events. 
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Figure 1: Top: Charged-particle yields for the different triggers normalized to the number of 
leading charged particles with pj > 0.4GeV/c in double-sided events, as a function 

of leading-track pj. The track-triggered distributions are normalized by the number of leading 
tracks in regions not affected by the rapid rise of the trigger efficiency near threshold. Bottom: 
Ratios of the leading-track pj distributions for the four different triggers. The stars indicate 
the ratio of the 12GeV/c over the minimum-bias samples, the circles the 20 over the 12GeV/c 
samples, and the squares the ratio of the 30 over the 20 GeV/c track-triggered spectra. 


Some events selected by the track triggers in Fig.j^are observed to result in a leading charged- 
particle Pj below the corresponding trigger threshold. This can happen if the t] of the track 
above threshold is outside the t] range considered in the analysis, and because the final track 
reconstruction—described in Section [0] —is more robust and selective than the fast-tracking al¬ 
gorithm implemented in the HLT. When the HLT selects an event based on a misreconstructed 
track, it is often the case that the track is not found in the final reconstruction. To determine 
the inclusive particle spectrum, events are first uniquely classified into leading-track pj classes 
in the pseudorapidity range in which the spectrum is being measured. The spectra are con¬ 
structed by taking events from the minimum-bias, 12 GeV/c track, 20 GeV/c track, and 30 GeV/c 
track trigger, respectively, for each bin. A 4% systematic uncertainty on the possible trigger-bias 
effect is estimated from MG simulations. This procedure was verified in a data-driven way by 
constructing a charged-particle spectrum from an alternative combination of event samples 
triggered by reconstructed jets. Both final spectra, triggered by tracks and jets, are found to be 
consistent within the associated systematic uncertainty. 

2.3 Track reconstruction 

The Pj distribution in this analysis corresponds to that of primary charged particles, defined as 
all charged particles with a mean proper lifetime greater than 1 cm/c, including the products of 
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strong and electromagnetic decays, but excluding particles originating from secondary interac¬ 
tions in the detector material. Weak-decay products are considered primary charged particles 
only if they are the daughters of a particle with a mean proper lifetime of less than 1 cm/c, 
produced in the collision. 

Charged particles are reconstructed using the standard CMS combinatorial track finder Il35l . 
The proportion of misreconstructed tracks in the sample is reduced by applying an optimized 
set of standard tracking-quality selections, as described in Ref. JSS]. A reconstructed track is 
considered as a primary charged-particle candidate if the statistical significance of the observed 
distance of closest approach between the track and the reconstructed collision vertex is less 
than three standard deviations, under the hypothesis that the track originated from this vertex. 
In case an event has multiple reconstructed collision vertices but is not rejected by the pileup 
veto, the track distance is evaluated relative to the best reconstructed collision vertex, defined 
as the one associated with the largest number of tracks, or the one with the lowest if multiple 
vertices have the same number of associated tracks. To remove tracks with poor momentum 
reconstruction, the relative uncertainty of the momentum measurement cr(px) / Px is required to 
be less than 10%. Only tracks that fall in the kinematic range of |piab| < 2.4 and px > 0.4GeV/c 
are selected for analysis to ensure high tracking efficiency (70-90%) and low misreconstruction 
rates (<2%). 

The yields of charged particles in each px arid p bin are weighted by a factor that accounts for 
the geometrical acceptance of the detector, the efficiency of the reconstruction algorithm, the 
fraction of tracks corresponding to a non-primary charged particle, the fraction of misrecon¬ 
structed tracks that do not correspond to any charged particle, and the fraction of multiply- 
reconstructed tracks, which belong to the same charged particle. 

The various correction terms are estimated using simulated minimum-bias pPb events from 
the HIJING event generator. To reduce the statistical uncertainty in the correction factors at 
high Px, samples of HijiNG events are also mixed with pp dijet events from the pythia MC 
generator II36H (version 6.423, tune D6T with CTEQ6L1 PDF for 2.76 TeV, tune Z2 for 7TeV Il37l l. 

The efficiency of the charged-particle reconstruction as well as the misreconstruction rates are 
also evaluated using pPb events simulated with EPOS. Differences between the two MC mod¬ 
els are mostly dominated by the fraction of charged particles consisting of strange and multi- 
strange baryons that are too short-lived to be reconstructed unless they are produced at very 
high Px. Such differences in particle species composition, which are largest for particles with 
3 < Px < 14GeV/c, are propagated as a systematic uncertainty in the measured spectra. Below 
this Px range, the strange baryons are only a small fraction of the inclusive charged particles in 
either model, and the difference in reconstruction efficiency between particle species has less 
impact at even larger px, as high-px multi-strange baryons can be directly tracked with high 
efficiency. Additional checks were performed by changing cutoffs imposed during track selec¬ 
tion and in the determination of the corresponding MG-based corrections. The corresponding 
variations in the corrected yields amount to 1.2-4.0% depending on the px region under con¬ 
sideration, and are included in the systematic uncertainty. 

Finite bin-widths and finite transverse momentum resolution can deform a steeply falling px 
spectrum. The data are corrected for the finite bin-width effect as they will be compared to a 
pp reference spectrum obtained by interpolation. The binning corrections are derived by fitting 
the measured distribution and using the resulting fit function as a probability distribution to 
generate entries in a histogram with the same px binning as used in the measurement. The cor¬ 
rection factors are then obtained from the ratio of entries in the bins of the histogram to the fit 
function evaluated at the centers of the bins. This correction amounts to 0-12%, depending on 
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Pj. a similar method is used to evaluate the "smearing" effect of the finite px resolution on the 
birmed distributions. It is found that the momentum measurement, which has a resolution of 
cr{pj)/pj « 1.5% near a px of lOOGeV/c, is sufficiently precise to only have a negligible effect 
on the measured spectra and therefore no correction factor is applied. To account for possi¬ 
ble incorrect determination of the momentum resolution from the simulation, the effects were 
again evaluated after increasing the value of cr{pj)/pj by an additional 0.01, which produces 
a maximal distortion in the spectrum at a given px of less than 1%. 

2.4 Proton-proton reference spectrum 

The pPb collisions occur at a center-of-mass energy of 5.02 TeV per nucleon pair. At this colli¬ 
sion energy, no proton-proton collisions have been provided by particle accelerators yet. The 
pp results closest in center-of-mass energy (y^) and with similar acceptance are those mea¬ 
sured at 2.76 and 7TeV by the CMS experiment Il20l l38ll . The determination of the nuclear 
modification factor Rab resides in an interpolated reference spectrum to be constructed from 
data at higher and lower energies. We follow the direct interpolation method developed in 
Ref. 1381 using measured px spectra from inelastic collisions with \rj\ < 1.0 at ^ys = 0.63, 1.8, 
and 1.96TeV collision energies from CDF 113^140^ , and 0.9, 2.76, and 7TeV collision energies 
from CMS l20H38l. This interpolation can be performed either as a function of px or as a func¬ 
tion of xj = 2pjc/^^/s. 

Since the px or xj values of the input data points are often different for each measurement 
performed at the various collision energies, each spectrum must first be fitted as a function 
of Px or xj. An interpolation is performed by fitting each of the spectra to a power-law de¬ 
pendence, and the resulting residuals to first- or third-order splines. The fitted spectra are 
then interpolated to determine the value of the reference spectrum at = 5.02 TeV using a 
second-order polynomial in the plane of the log-log invariant production vs. ^/s, as shown in 
Fig II For the px-based direct interpolation, data from only two of the six spectra are avail¬ 
able at Px > 30GeV/c, which implies that the px-based direct interpolation is well constrained 
only at low px- On the other hand, the vx-based interpolation is well constrained at high px 
for y/s = 5.02 TeV, so it is natural to combine the reference distributions from these two direct 
interpolation methods. 

The final pp reference spectrum is obtained by combining the px- and vx-based reference spec¬ 
tra as follows. The px-based reference is chosen for px below 12.5GeV/c, and the Vx-based 
result above 13.5 GeV/c; between these two px values a linear weighting is implemented for the 
two references. The systematic uncertainty in the pp reference spectrum is determined through 
changing both the specific method of interpolation, as well as the underlying pp reference data 
within their statistical and systematic uncertainties. The systematic uncertainty is dominated 
by the interpolation method, and is determined by comparing the combined px- and Vx-based 
reference spectra to the reference spectra obtained solely from the px or xj distributions, and 
also from a reference spectrum determined by a "relative placement" method. In the latter, the 
reference spectrum is obtained by computing where the 5.02 TeV spectrum is situated with re¬ 
spect to the 2.76 and 7 TeV spectra in pythia, and applying the computed placement factors to 
the measured 2.76 and 7 TeV spectra. The placement factors are determined by taking the value 
of the 5.02 TeV pythia spectrum, subtracting the value of the 2.76 TeV spectrum, and divid¬ 
ing by the difference between the 7 and the 2.76 TeV spectra. This process is then reversed by 
using the computed placement factors from pythia, and replacing the 2.76 and 7 TeV PYTHIA 
spectra with the measured ones to determine the interpolated 5.02 TeV spectrum. Addition¬ 
ally, the NLO-based center-of-mass energy rescaling proposed in Ref. BTI is found to yield 
results consistent within the uncertainties of the other employed methods. The uncertainty 
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Figure 2: Examples of interpolations between measured charged-particle differential cross sec¬ 
tions at different for pj values of 3 and 15 GeV/c (top left and right), xj values of 0.01 and 
0.02 (bottom left), and Xj values of 0.03 and 0.04 (bottom right). These Xj values correspond to 
Pj ~ 25, 50, 75, and lOOGeV/c at ^/s = 5.02 TeV. The second-order pol}momial fits, performed 
in the plane of the log-log invariant production vs. ^/s, are shown by the solid lines. The open 
squares and circles, and the filled crosses represent interpolated cross section values at 5.02 TeV 
using different methods: px'based interpolation, Vj-based interpolation, and relative place¬ 
ment, respectively. The error bars on the interpolated points represent the uncertainties in the 
fit. 


in the pp reference distribution due to the interpolation method is estimated to amount to 
10%, which captures the overall point-to-point variations in all of the interpolation and scaling 
methods employed. The contribution from the uncertainties in the underlying pp input data 
corresponds to 6%. These numbers are added in quadrature, resulting in the 12% uncertainty 
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quoted for the \/s = 5.02 TeV interpolated pp reference spectrum. 


3 Systematic uncertainties 

A summary of all the contributions to systematic uncertainties in the px spectra, Rppb, and Yasym 
are given in Table The asterisk symbol is introduced to denote that an interpolated, rather 
than measured, pp reference spectrum is used to construct the nuclear modification factor. 
Aside from the uncertainty from the spectra relative normalization and average nuclear thick¬ 
ness, all uncertainties are determined by taking the approximate maximum deviation from the 
central value found for the given source. For the particle species uncertainty, an asymmetric 
uncertainty band is quoted because the maximum deviation above the central value of the mea¬ 
surement is much larger than the maximum deviation below. For the purpose of determining 
the significance of observed features in the measurement, the uncertainties are conservatively 
treated as following a Gaussian distribution with a root mean square given by the value of the 
uncertainty as determined above. 

The degree of correlation among different uncertainties is described next. For the spectra and 
i?*Pb measurements, the uncertainty in the efficiency of the single-track trigger and offline 
requirements in selecting DS events is largely a normalization uncertainty, although it also 
slightly affects the shape of the spectrum for pi < 3 GeV/c. The uncertainty from the contribu¬ 
tion of the various particle species to the unidentified spectrum has the most significant effect 
in the region 3 < pj < 14GeV/c and can impact the shape of the spectrum in a smooth fash¬ 
ion. At high px, this effect is less prominent because, due to time dilation, unstable particles 
have a higher probability of traversing the irmer tracker before decaying and therefore a higher 
probability of being reconstructed. Therefore, from this uncertainty the lower bound on the 
pPb spectra measurement at higher px is 2.5% below the central value, which corresponds to 
no unstable particles being produced. Uncertainty in track misreconstruction can also affect 
the shape of the measured spectrum, as the misreconstructed fraction of high-px particles is 
sensitive to large occupancy in the silicon tracker within the cones of high-energy jets. The un¬ 
certainty in tracking selection can also affect the shape of the spectrum by raising or lowering 
the measured values at high px, without changing the low-px values, as high-px tracks are more 
sensitive to possible mismodeling of detector alignment than low-px tracks. The uncertainty 
in the relative normalization of spectra is computed from the normalization factors involved in 
the combination of the px distributions from different triggers. This uncertainty only applies 
for selected px regions, and may raise or lower the spectrum above px = 14 GeV/c by a constant 
factor of 1% relative to the lower-px part of the spectrum. The uncertainty from potential biases 
of the method used to combine triggers can also affect the shape of the spectrum above px = 
14GeV/c. 

For the measurement, the uncertainty in the average nuclear thickness function 13 can 
influence the curve by a constant multiplicative factor. The uncertainty from the pp inter¬ 
polation is strongly correlated among points close together in px, while some partial correlation 
remains throughout the whole px region, even for very different px values. 

For the forward-backward asymmetry measurements, most of these uncertainties cancel in 
part or in full when the ratio of the spectra is taken. The remaining uncertainty in the detector 
acceptance and tracking efficiency can change the shape of the forward-backward asymmetry, 
particularly at high px. 
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4 Results 


Table 1: Systematic uncertainties in the measurement of charged-particle spectra, and 

^asym- The ranges quoted refer to the variations of the uncertainties as a function of pj. Values 
in parentheses denote the negative part of the asymmetric uncertainty where applicable. The 
total uncertainties of the measured pPb and the interpolated pp spectra, as a function of pj, are 
shown in the lower panel of Fig.[^ _ 

Source Uncertainty [%] 

Trigger efficiency 1.0 

Momentum resolution 1.0 


Particle species composition 

1-10.0 (0.5-5) 

Track misreconstruction rate 

1.0 

Track selection 

1.2-4.0 

Spectra relative normalization 

O.O-l.O 

Trigger bias 

0.0-4.0 

Total (spectra) 

2.2-10.9 

pp interpolation 

12.0 

Total (R;p,) 

12.2-16.2 

(TpPb) average nuclear thickness 

4.8 

Total (Yasym 0.3 < //cm < 0.8) 

2.0-3.0 

Total (Yasym 0.8 < |//cm| < 1-3) 

2.0-5.0 

Total (Yasym 1-3 < |//cm| < 1-8) 

2.0-5.0 


4 Results 

The measured charged-particle yields in double-sided pPb collisions at = 5.02 TeV are 

plotted in Fig. I^for the IpcmI < 1.0,0.3 < ±t]cM < 0.8,0.8 < ±?/cm < 1.3, and 1.3 < ±ricM < 1-8 
pseudorapidity ranges. Positive (negative) pseudorapidity values correspond to the proton 
(lead) beam direction. To improve the visibility of the results, the spectra at different pseu¬ 
dorapidities have been scaled up and down by multiple factors of 4 relative to the data for 
I 7cm I < 1- The relative uncertainties for the pPb and the pp spectra are given in the bottom 
panel. 

The measurement of the charged-particle nuclear modification factor of Eq. Q requires a rescal¬ 
ing of the pp cross section by the average nuclear thickness function in minimum-bias pPb 
collisions. This factor amounts to (Tppb) = (0.0983 ± 0.0044) mb^^ for inelastic pPb colli¬ 
sions and is obtained from a Glauber MC simulation EM, where the Pb nucleus is de¬ 
scribed using a Woods-Saxon distribution with nuclear radius 6.62 ± 0.13 fm and skin depth 
of 0.546 ± 0.055 fm Il3lll3l. As double-sided events correspond to 94-97% of inelastic collisions 
based on HIJING and EPOS MG computations 1311, the value of (TpPb) would be about 5% 
higher for double-sided events. 

The charged-particle Rppb at mid-rapidity (|7 cm| < 1) is plotted in Fig. |4 as a function of pj. 
The shaded band at unity and pj ~ 0.6 represents the uncertainty in the Glauber calculation of 
(Tppb). The smaller uncertainty band around the measured values shows the fully correlated 
uncertainties from the following sources: spectra relative normalization, track selection, and 
trigger efficiency. The total systematic uncertainties are shown by the larger band around the 
measured values (Table [^. The nuclear modification factor shows a steady rise to unity at 
Pt ~ 2 GeV/c, then remains constant at unity up to approximately 20 GeV/c, and rises again at 
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Figure 3: Top: Measured charged-particle transverse momentum spectra in pPb collisions at 
= 5.02 TeV for: \rjcu\ < 1-0, 0.3 < ±/ycM < 0.8, 0.8 < ±?/cm < 1-3, and 1.3 < ±?/cm < 1-8, 
and the interpolated pp reference spectrum in | //cm | < 1, normalized to the number of double¬ 
sided events. Positive pseudorapidity values correspond to the proton beam direction. The 
spectra have been scaled by the quoted factors to provide better visibility. Bottom: Systematic 
uncertainties in the measured pPb and interpolated pp spectra, as a function of pj (see text). 


higher pj, reaching a maximum value around 1.3-1.4 above 40 GeV/c. 

The fact that the nuclear modification factor is below unity for pj < 2 GeV/c is anticipated 
since particle production in this region is dominated by softer scattering processes, that are not 
expected to scale with the nuclear thickness function. In the intermediate pj range (2-5 GeV/c), 
no significant deviation from unity is found in the Rpp^ ratio. There are several prior mea¬ 
surements that suggest an interplay of multiple effects in this pj region. At lower collision 
energies, an enhancement ("Gronin effect" II44I ') has been observed IITStiTSi that is larger for 
baryons than for mesons, and is stronger in the more central collisions. This enhancement has 
been attributed to a combination of initial-state multiple scattering effects, causing momen¬ 
tum broadening, and hadronization through parton recombination (a final-state effect) [45l 
invoked to accommodate baryon/meson differences. Recent results from pPb collisions at 
~ 5.02TeV llMll46ti4^ and from dAu collisions at = 200GeV IISOllMlI suggest that 
collective effects may also play a role in the intermediate-pj region. Most theoretical models do 
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4 Results 


not predict a Cronin enhancement in this pj range at LHC energies as the effect of initial-state 
multiple scattering is compensated by nPDF shadowing Il52l . 

In Fig.|^ the CMS measurement is compared to the result of an NLO pQCD calculation Il5^ for 
charged particles produced at mid-rapidity The calculation uses the CTEQIO 15^ free-proton 
FDF, the EPS09 nPDF |^, and the fDSS fragmentation functions 1551 . The observed rise of the 
nuclear modification factor up to Rppt ~ 1.3-1.4 at high pj is stronger than expected theoret¬ 
ically None of the available theoretical models Il5^ predict enhancements beyond Rppb ~ 1.1 
at high Pj. In particular, although the pj range corresponds to parton momentum fractions 
0.02 < ^ < 0.2, which coincides with the region where parton antishadowmg effects are ex¬ 
pected IflOl . none of the nPDFs obtained from global fits to nuclear data predict enhancements 
beyond 10% at the large virtualities (Q^ p^ 500-10 000 GeV/c^) of relevance here. 

An estimate of the significance of this observed rise above unity for 40 < px < 120 GeV/c is 
determined by interpreting all uncertainties as following a multivariate normal distribution 
where the components are the six pj bins in the kinematic region of interest. The variance of 
each component is given by the sum of the statistical and systematic uncertainties in quadra¬ 
ture. For the case of the asymmetric particle species uncertainty, the smaller negative value is 
used as the data are uniformly larger than the expected values of the hypothesis to be tested. 
Given that the uncertainties of the reference spectrum are derived from applying different in¬ 
terpolation procedures and propagating the uncertainties from previous measurements from 
multiple experiments, it is not possible to unambiguously determine how all systematic un¬ 
certainties are correlated between measurements in each pj bin. Therefore, a pair of estimates 
of the possible significance is given. In one case, only the systematic uncertainties from the 
relative normalization of the spectra, track selection, trigger efficiency, nuclear thickness func¬ 
tion, and NLO pQGD calculation are treated as fully correlated, while others are treated as 
uncorrelated. In the other case, all systematic uncertainties are treated as fully correlated. Both 
the hypothesis that Rpp^ is unity and the hypothesis that Rppj, is given by the NLO pQGD 
calculation are tested. For the case in which some uncertainties are treated as uncorrelated, a 
log-likelihood ratio test is performed using an alternative hypothesis that Rppb is given by ei¬ 
ther unity or the NLO prediction, scaled by a constant, pj-independent, factor. The hypothesis 
that Rppb is unity for 40 < px < 120 GeV/c is rejected with a p value of 0.006%, and the hypoth¬ 
esis that Rppb is given by the NLO pQGD calculation for 40 < px < 120 GeV/c is rejected with 
a p value of 0.2%. For the case in which all uncertainties are fully correlated, the log-likelihood 
ratio test cannot be used, as the covariance matrix becomes nearly singular and the maximum 
likelihood estimation fails. Instead, a two-tailed univariate test is performed using the single 
measurement for 61 < px < 74GeV/c. From this test, the hypothesis that Rppb is unity for 
61 < Px < 74GeV/c is rejected with a p value of 0.4%, and the hypothesis that Rppb is given by 
the NLO pQGD calculation for 61 < px < 74 GeV/c is rejected with a p value of 2%. 

Fig. 1^ also shows the measurement from the ALIGE experiment |i22l , which is performed in a 
narrower pseudorapidity range than the GMS one, and uses a different method (NLO scaling) 
to obtain the pp reference spectrum based on ALIGE pp data measured at ^/s = 7 TeV. The 
difference in the GMS and ALIGE Rppi, results stems primarily from differences in the charged- 
hadron spectra measured in pp collisions at = 7TeV Il38ll56ll . 

Figure shows the forward-backward yield asymmetry, Yasym (Eq. |^, as a function of px for 
0.3 < ^cm| < 0-8, 0.8 < IpcmI < 1-3, and 1.3 < |Pcm| < 1-8. In all three rj ranges, the 
value of Yasym rises from px ~ 0.4 to about 3GeV/c, then falls to unity at a px of lOGeV/c, 
and remains constant at unity up to the highest px values. At the lowest px value, Yasym is 
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Figure 4: Measured nuclear modification factor as a function of px for charged particles pro¬ 
duced in IpcmI < 1- The shaded band at unity and pj « 0.6 represents the uncertainty in the 
Glauber calculation of (TpPb). The smaller uncertainty band around the data points shows the 
uncertainty from effects (combining spectra, track selection, and trigger efficiency) that are fully 
correlated in specific px regions. The total systematic uncertainties, dominated by uncertainty 
in the pp interpolation, are shown by the larger band (see Table [^. 
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consistent with unity for 0.3 < |Pcm| < 0.8, but is above unity in the larger pseudorapidity 
regions. For px < lOGeV/c, the Yasym is larger than unity as has been predicted by models 
including nuclear shadowing [52|. A theoretical NLO pQGD computation of Yasym at high 
Px l53l , using GTEQ6 USTl free-proton PDFs, EPS09 nPDFs |4i|, and Kretzer parton-to-hadron 
fragmentation functions Il58ll , is also shown in Fig.The theoretical predictions are consistent 
with these data. 

To determine if the Rppb and Yasym results can be consistently interpreted in terms of nPDF 
modifications, an MG study using the pythia (Z2 tune) event generator was performed to 
correlate each high-px hadron to the fractional momentum, x, of the initial-state parton from 
the Pb nucleus that participated in the hard-scattering process producing the final hadron. In 
all pseudorapidity intervals studied here, most of the hadrons with px > 20 GeV/c, i.e., in the 
range where the Rpp^ exceeds unity in Fig. |4j come from the x region that is associated with 
antishadowing in the nPDF distributions. Although the mean of the x distribution increases 
with rjcM, for hadrons with px above 20GeV/c it remains in the range 0.02 < v < 0.2. Thus, 
similar antishadowing effects are expected in the positive and negative rjcu regions resulting in 
a Yasym close to unity. At low px, corresponding to v < 0.02, a larger hadron yield is observed 
in the direction of the Pb beam. This is qualitatively consistent with expectations of gluon 
shadowing 1521 . 

An enhancement in Rppj, at high px can possibly arise if the quark-jet fraction is larger in pPb 
than in pp collisions. Since the charged-particle products of quark fragmentation more often 
have higher relative px than those produced by gluon fragmentation, that could lead to an en- 
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5 Summary 



Figure 5: Charged-particle nuclear modification factors measured by CMS in |?/cm| < 1 (filled 
circles), and by ALICE in \rjcu\ < 0.3 (open squares), are compared to the NLO pQCD pre¬ 
diction of Ref. I53l . The theoretical uncertainty is based on the EPS09 error sets. For the CMS 
measurement, the shaded band at unity and px ~ 0.6 represents the uncertainty in the Glauber 
calculation of (Tppb), the smaller uncertainty band around the data points shows the fully cor¬ 
related uncertainties and the total systematic uncertainty is shown by the larger band (see Ta¬ 
ble]^. For the ALICE measurement, the total systematic uncertainties, excluding the normal¬ 
ization uncertainty of 6%, are shown with open boxes. 


hancement in the charged-particle production at high px beyond NLO expectations, without 
a corresponding increase in the jet Rppb Il25l[26l . We note that the gluon-to-hadron fragmen¬ 
tation functions are not well constrained in pp collisions at LHC energies 1271 . although such 
uncertainties should largely cancel in ratios of cross sections. 


5 Summary 

Charged-particle spectra have been measured in pPb collisions at = 5.02 TeV in the trans¬ 
verse momentum range of 0.4 < px < 120GeV/c for pseudorapidities up to |Pcm| = 1-8. The 
forward-backward yield asymmetry has been measured as a function of px for three bins in tjcu- 
At Px < 10 GeV/c, the charged-particle production is enhanced in the direction of the Pb beam, 
in qualitative agreement with nuclear shadowing expectations. The nuclear modification factor 
at mid-rapidity, relative to a reference spectrum interpolated from pp measurements at lower 
and higher collision energies, rises above unity at high px reaching an Rpp^, value of 1.3-1.4 at 
Ft ^ 40 GeV/c. The observed enhancement is larger than expected from NLO pQGD predic¬ 
tions that include antishadowing effects in the nuclear PDFs in this kinematic range. Future 
direct measurement of the spectra of jets and charged particles in pp collisions at a center-of- 
mass energy of 5.02 TeV is necessary to better constrain the fragmentation functions and also 
to reduce the dominant systematic uncertainties in the charged-particle nuclear modification 




















15 



[GeV/c] 


Figure 6: Charged-particle forward-backward yield asymmetry as a function of px for 0.3 < 
IPcmI < 0.8 (top), 0.8 < IpcmI < 1-3 (middle), and 1.3 < |Pcm| < 1-8 (bottom). The asymmetry 
is computed as the charged-particle yields in the direction of the Pb beam divided by those of 
the proton beam. The solid curves are NLO pQCD theoretical calculations including nPDFs 
modifications 1531 . The theoretical uncertainty is based on the EPS09 error sets. 


factor. 
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